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(57) Abstract: The present invention provides a new method and device for disease detection, more particularly cancer detection, 
from the analysis of diffuse reflectance spectra measured in vivo during endoscopic imaging. The measured diffuse reflectance spec- 
tra are analyzed using a specially developed light-transport model and numerical method to derive quantitative parameters related to 
tissue physiology and morphology. The method also corrects the effects of the specular reflection and the varying distance between 
endoscope tip and tissue surface on the clinical reflectance measurements. The model allows us to obtain the absorption coefficient 
(jLia) and further to derive the tissue micro-vascular blood volume fraction and the tissue blood oxygen saturation parameters. It also 
allows us to obtain the scattering coefficients (jus and g) and further to derive the tissue micro -particles volume fraction and size 
distribution parameters. 
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Method And Apparatus For Measuring Cancerous Changes From Reflectance Spectral 
Measurements Obtained During Endoscopic Imaging 

5 

Inventors: Yasser Sherif Fawzy, a citizen of Egypt 
Haishan Zang, a citizen of Canada 

1 0 BACKGROUND OF THE INVENTION 

The present invention relates to the field of optical spectroscopy and more particularly 
to the method for obtaining information about tissue physiology and morphology using 
diffuse reflectance spectroscopy. The purpose of the invention is to develop a non-invasive 
optical method for cancer detection. 

15 Lung cancer is the leading cause of cancer death in North America, and it has the 

second most common cancer incidence among both men and women. Medical research 
indicates that cancer can be treated more effectively when is detected early, when lesions are 
smaller or when tissue is in a precancerous stage. Unfortunately, conventional lung 
endoscopy (bronchoscopy) based on white light reflectance (WLR) imaging, which is used 

20 typically to detect the cancer lesions in the central airways of the lung, can only detect about 
25 percent of the lung cancers. Most of these lesions are in the late stage when cancer has 
progressed and is fatal. This detection rate has created the need for a detection or imaging 
modality to accompany WLR imaging and achieve better diagnostic performance for cancer 
detection. 

25 A number of research groups have investigated the use of tissue autofluorescence to 

improve the detection sensitivity of cancerous lesions. Just as certain morphological changes 
in tissue may be associated with disease, chemical changes may also be exploited for disease 
detection especially for early detection of disease. When tissue is illuminated (or excited) 
with specific wavelengths of ultraviolet (UV) or visible light, biological molecules 

30 (fluorophores) will absorb the energy and emit it as fluorescent light at longer wavelengths 
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(green/red wavelength region). These wavelengths of light are selected based on their ability 
to stimulate certain chemicals in tissue that are associated with disease or disease processes. 
Images or spectra from these emissions (fluorescence) may be captured for observation 
and/or analysis. Diseased tissue has considerably different fluorescence signals than healthy 
5 tissue so the spectra of fluorescence emissions can be used as a diagnostic tool. 

In United States Published Patent Application No. 2004/245350 to Zeng, entitled 
"Methods and Apparatus for Fluorescence Imaging using Multiple Excitation-Emission Pairs 
and Simultaneous Multi-Channel Image Detection", the inventor reports use of a second 
independent fluorescence signal in the red/NIR wavelength region. The diseased tissue such 

10 as cancerous or pre-cancerous tissue illuminated with the red/NIR light, unlike the tissue 
properties discussed above, emits fluorescence, providing intensities that are higher for 
diseased tissue than for normal tissue. These properties may be exploited to improve image 
normalization and diagnostic utility of images. 

Although fluorescence imaging provides increased sensitivity to diseases such as 

15 cancer, there are also some trade offs. A commercial fluorescence imaging system has 
achieved sensitivity of 67 percent for lung cancer detection. However, such increase in 
detection sensitivity was at the cost of the decreased detection specificity, which was reduced 
to 66 percent compared to 90 percent for WLR imaging alone. The result was increased 
medical costs related to the enlarged number of biopsies caused by the increased number of 

20 false positives. 

In order to provide more accurate diagnosis of cancerous tissue, a more convenient 
approach has been to perform additional non-invasive and real-time cancer diagnosis that 
would increase detection specificity, reduce medical cost, and help doctors during surgery to 
define cancerous region of the tissue. There are few known methods of non-invasive cancer 
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diagnosis, such as reflectance spectroscopy and fluorescence spectroscopy, both of which are 
based on detection of biochemical and morphological variation of the diseased tissue. 

Biological tissue is a turbid medium, which absorbs and scatters incident light. When 
light impinges on tissue, it is typically multiple elastically scattered but at the same time 
5 absorption and fluorescence can occur, too. Further scattering and absorption can occur 
before light exits the tissue surface containing compositional and structural information of the 
tissue. This information can be used for detection of pre-cancers and early cancers that are 
accompanied by local metabolic and architectural changes at the cellular and subcellular 
level, for example, changes in the nuclear-to-cytoplasm ratio of cells and changes in 

10 chromatin texture. These changes affect the elastic scattering properties of tissue. 

Reflectance spectroscopy is an analysis of a light reflected from tissue. Tissue 
reflectance spectroscopy can be use to derive information about tissue chromophores 
(molecules that absorbs light strongly), e.g. hemoglobin. The ratio of oxyhemoglobin and 
deoxy-hemoglobin can be inferred and used to determine tissue oxygenation status, which is 

15 very useful for cancer detection and prognosis analysis. It can also be used to derive 
information about scatterers in the tissue, such as the size distribution of cell nucleus and 
average cell density. In many cases quantification of chromophore concentration is desired, 
and this requires the ability to separate the effects of absorption from those of scattering. 

Fluorescence spectroscopy is the analysis of fluorescence emission from tissue. 

20 Native tissue fluorophores (molecules that emit fluorescence when excited by appropriate 
wavelengths of light) include tyrosine, tryptophan, collagen, elastin, flavins, porphyrins and 
nicotinamide adenine dinucleotide (NAD). Tissue fluorescence is very sensitive to chemical 
composition and chemical environment changes associated with disease transformation. 
Exogenous or exogenously-induced chromophores that have been shown to accumulate 

25 preferentially in the diseased areas can also be used. 
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Another type of spectroscopic technique that has been used to examine tissues 
involved the use of Raman spectroscopy. Raman spectra convey specific information about 
the vibrational,, stretching, and breaking bond energies of the illuminated sample. Raman 
spectroscopy probes molecular vibrations and gives very specific, fingerprint-like spectral 
5 features and has high accuracy for differentiation of malignant tissues from benign tissues. 
Raman spectroscopy can also be used to identify the structural and compositional differences 
on proteins and genetic materials between malignant tissues, their pre-cursers, and normal 
tissues. The development of an in vivo tissue Raman probe, however, is technically 
challenging due to the weak Raman signal of tissue, interference from tissue fluorescence, 

10 and spectral contamination caused by the background Raman and fluorescence signals 
generated in the fiber itself. 

Another non-invasive imaging technology is optical coherence tomography (OCT). It 
is based on the principle of low-coherence interferometry where distance information 
concerning various tissue microstructures is extracted from time delays of reflected signals. 

15 OCT can perform high-resolution "optical biopsies" of tissue microstructure in situ and in 
real-time. However, the spatial resolution of commercial OCT systems still cannot meet the 
clinical requirements for accurate in vivo endoscopy diagnosis. 

Other than these methods, there are on-going research studies in the area of non- 
invasive cancer diagnosis based on morphological alteration of cell structure for cancer cells. 

20 One of the most prominent features used by pathologist to diagnose tissue as being cancerous 
is the presence of enlarged and crowded nuclei. Since the nuclei of cancerous cells are 
significantly larger than nuclei of normal cells for many cancer types, the target of these 
research studies is to estimate the average size of scatterers such as nuclei, mitochondria, and 
other organelles of cells, non-invasively through an optical system. 
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When a beam of light reaches the tissue under investigation, part of it will be 
specularly reflected by the surface, while the rest is refracted and transmitted into the tissue. 
The light transmitted into the tissue will be scattered and absorbed. After multiple scattering, 
some of the transmitted light will return to the tissue surface and appear for detection. Light 
5 scattering by biological tissues are caused by refractive index variations inside the tissue at 
the boundaries of various microstructures such as cell nucleus and collagen bundles. Thus, 
tissue scattering property changes with variations in a tissue's micro structure properties and 
morphology, which are often accompanied with tissue pathological changes. For example, 
when normal tissue becomes cancerous, the nucleus size of the cells and the epithelial layer 

10 thickness increase as does the total volume occupied by the cells (micro - s catterer s) . Such 
changes in the tissue microstructure and morphology have been found to cause intrinsic 
differences in the light-scattering properties of the normal and cancerous lesions. 

In particular, two measurement approaches could be identified in literature for 
obtaining quantitative differences on scattering properties of normal and cancerous lesions 

15 using reflectance spectral measurements. One approach is to measure the single light- 
scattering spectra (LSS) originated from the superficial tissue layers, and to extract 
quantitative information about the scattering structures at the cellular and sub-cellular levels. 
The LSS technique examines variations in the elastic scattering properties of cell organelles 
to infer their sizes and other dimensional information. In order to measure cellular features in 

20 tissues and other cellular structures, it is necessary to distinguish the weak, singly scattered 
light from diffuse light, which has been multiple scattered and no longer carries easily 
accessible information about the scattering objects. Consequently, concentration of the 
scatterers in suspension must be low so that information obtained from only angular 
distribution of single scattered photons can be analyzed. Analysis of the singly-backscattered 
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light spectrum using light-scattering theory provides information about the size and number 

density of cell nuclei without tissue removal. 

Nevertheless, these LSS measurements are limited since LSS does not allow obtaining 

quantitative information about the absorption properties of the tissue such as chromophore 

5 concentration. In addition, LSS measurements are difficult if not impossible to perform 

during endoscopy applications. 

The other approach is to obtain quantitative information about tissue morphology 

(tissue scattering properties) from the diffuse reflectance spectra (DRS). Diffuse reflectance 

relies upon the projection of a light beam into the sample where the light is reflected, 

10 scattered, and transmitted through the sample material. The back-reflected, diffusely 

scattered light (some of which is absorbed by the sample) is then collected by the accessory 

and directed to the detector optics. Only the part of the beam that is scattered within a sample 

and returned to the surface is considered to be diffuse reflection. 

Diffuse reflectance measurements are simpler to implement and allow obtaining 

15 quantitative information about the absorption properties as well as the scattering properties. 

However, in most studies quantitative information obtained from DRS measurements was 

limited to the estimation of the average bulk tissue optical properties (reduced scattering and 

absorption coefficient) rather than obtaining quantitative information related directly to tissue 

microstructure and morphology. This limitation is mainly due to the complex nature of light 

20 propagation (multiple scattering) in tissue with such microstructures and morphology. 

Therefore, it is difficult to characterize the scattering properties at cellular levels from DRS. 

Also, the back-reflected light can be considered as derived from two categories, 

diffuse reflectance and specular reflectance. The specular reflectance is the light that does 

not propagate into the sample, but rather reflects from the front surface of the tissue. This 

25 component contains information about the tissue at the surface. The diffuse component is 
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generally considered more useful for tissue qualification and quantification than is the 

specular component. 

Various approaches, such as using a contact fiber-optic probe, collecting returning 
radiation over a small collection angle, or using a specular control device, have been 
5 proposed to emphasize the diffuse component relative to the specular component. For some 
tissues, for example, skin, it is relatively easy to obtain such spectra by simply touching the 
lesion with an appropriate optical fiber bundle that is coupled to a spectrometer. However, 
for internal organs such as the lung, such set-up would not be practical because of the 
interferences of the instrument-channel-based fiber probe with biopsy or other therapeutic 
10 tools. 

Few studies have investigated the potential of diffuse reflectance spectroscopy for 
detecting tissue cancerous changes. Intrinsic differences on optical properties between 
malignant and benign lesions/normal tissues were found and were related directly to the 
changes in tissue physiology and morphology that occurred during cancer transformation. 

15 Clinical spectroscopic measurements and analyses have been performed on various organ 
sites including the lung. In particular. Bard et ah have performed spectral measurements and 
analysis on the abnormal lesions that were identified during fluorescence bronchoscopy and 
they found significant changes on both absorption-related and scattering-related physiological 
and morphological properties when tissue became malignant. They have also evaluated the 

20 potential of such spectral measurements for improving lung cancer detection specificity. 
However, their measurements were still conducted using a fiber optic probe inserted through 
the endoscope instrument channel and been in contact with the tissue surface during the 
measurement. 

In principle, DRS as used in the clinical setting is performed in the following manner. 
25 An optical fiber probe, fiber-optic bundle, inserted through the biopsy channel of the 



7 



WO 2006/076810 PCT/CA2006/000080 

endoscope and, coupled to a spectrometer, is brought into contact with the tissue surface. The 
optical fiber probe consists of an illuminating fiber/fiber optical bundle, typically the central 
core and surrounding fibers/ fiber optical bundles for capturing the returning radiation. Light 
leaves the illuminating fiber and enters the tissue under investigation. After the processes of 
scattering and absorption, light that leaves the tissue is captured by the detecting fibers and 
directed into a spectrometer. The spectrum may than further analyzed to determine the 
characteristics of the tissue. 

Despite the fact that employment of the contact probe geometry gives more controlled 
diffuse reflectance measurements with less measurement artefacts, the limitation of this kind 
of measurement is that it is awkward and time consuming for in vivo endoscopic imaging of 
internal organs. 

There is therefore a need for an apparatus to obtain quantitative information about 
tissue physiological and morphological characteristics directly from diffuse reflectance 
measurements with the goal for endoscopy applications. Accordingly, the present invention 
uses a non-contact probe, which eliminates the need for a fiber probe through the instrument 
channel. Thus, the present invention overcomes the problems presented in the prior art and 
provides additional advantages over the prior art. 

BRIEF SUMMARY OF THE INVENTION 

The present invention in one embodiment is a method of obtaining information about 
tissue from diffuse reflectance spectra measured in vivo during endoscopic imaging by 
illuminating a tissue with a broadbeam radiation to produce returning radiation; measuring a 
diffuse reflectance spectra of the returning radiation with a non-contact probe; analyzing the 
diffuse reflectance spectra for a two-layer tissue model by one-dimensional light 
transportation modeling; extracting at least one optical property of the tissue from the 
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analyzed diffuse reflectance spectra; and deriving information about at least one of a 
physiology and a morphology of the tissue from the optical property. 

In another embodiment, the present invention is an apparatus for obtaining 
information about tissue from diffuse reflectance spectra, having means for illuminating a 
5 tissue with a broadbeam radiation to produce returning radiation; a non-contact probe to 
measure the returning radiation; means for measuring a diffuse reflectance spectra of the 
returning radiation; means for analyzing the diffuse reflectance spectra for a two-layer tissue 
model by one-dimensional light transportation modeling; means for extracting at least one 
optical property of the tissue from the analyzed reflectance spectra; and means for deriving 
10 information about at least one of a physiology and a morphology of the tissue from the 
optical property. 

In yet another embodiment, the present invention is an apparatus for obtaining 
information about tissue from diffuse reflectance spectra, having a non-contact probe; a light 
source producing a broadband interrogating radiation at a distal end of said probe to 
15 illuminate a tissue and to produce returning radiation; a detecting system coupled to capture 
said returning radiation; a processing unit coupled to said detecting system, said processing 
unit measuring a diffuse reflectance spectra of said returning radiation and classifying the 
tissue as one of benign and malignant based on said diffuse reflectance spectra. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

The organization and manner of the structure and operation of the invention, together 
with further objects and advantages thereof, may best be understood by reference to the 
following description, taken in connection with the accompanying drawings, wherein like 
25 reference numerals identify like elements in which: 
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FIG. 1 is an illustration of an endoscopy system of the preferred embodiment of the 
present invention; 

FIG. 1 A is a cross-sectional view of a probe used in the endoscopy system of FIG. 1; 
FIG. IB is a diagram of a spectral attachment used in the endoscopy system of FIG. 1; 
5 FIG. 2 is a diagram of the measuring geometry; 

FIG 2a is a graph showing a light fluence distribution ^ as a function of tissue depth z; 
FIG. 3 is a diagram of normalization procedure; 
FIG. 4 is a diagram of forward model; 

FIG. 5 is a flow diagram of the procedure for disease detection based on the 
10 information obtained from the analysis of the diffuse reflectance spectra; 

FIG. 6 contains graphs of (a) reflectance spectra measured and fitted (lines) from two 
normal tissue sites/benign lesions and two malignant lung lesions from the same patient; and 
(b) corrected true diffuse reflectance spectra derived from Eq. (3) using the correction 
parameters ao and bo obtained with the developed model; 
15 FIG. 7 is a set of graphs showing absorption and scattering related parameters 

obtained from fitting the two benign and two malignant spectra: (a) Blood volume fraction 
(p) 5 (b) Tissue oxygen saturation parameter (a), (c) Scattering volume fraction in both layers 
(51 and 52), and (d) Size-distribution parameter describing the scatterer's size-distribution in 
both layers (pi and 02); 

20 FIG. 8 contains graphs of the average reflectance spectra fitted and corrected obtained 

from the 50 normal tissue/benign lesions versus that of the 50 malignant lesions; 

FIG. 9 contains scatter plots of the physiological and morphological parameters of the 
bronchial mucosa obtained from the reflectance spectral analysis: (a) Blood volume fraction p 

and (b) Oxygen saturation parameter a, (c) Scattering volume fraction and (d) Size- 
25 distribution parameter P 1 ; and 

10 
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FIG. 10 is a binary plot of (a) blood volume fraction versus the scattering volume, and 
(b) blood volume fraction versus oxygen saturation parameter. 



11 



WO 2006/076810 PCT/CA2006/000080 

DETAILED DESCRIPTION OF THE INVENTION 

While the invention may be susceptible to embodiments in different forms, there is 
shown in the drawings, and herein will be described in detail, specific embodiments with the 
understanding that the present disclosure is to be considered an exemplification of the 
5 principles of the invention, and is not intended to limit the invention to that as illustrated and 
described herein. 

The approach of the present invention is to obtain quantitative information about the 
absorption-related and/or scattering-related properties from the diffuse reflectance spectra 
(DRS) obtained during in vivo endoscopic imaging. Diffuse reflectance relies upon the 

10 projection of a broadband light beam into the sample where the light is absorbed, reflected, 
scattered, and transmitted or back-reflected through the sample material. The back-reflected 
(back-scattered) light is then collected by the accessory (e.g. an optical fiber) and directed to 
the detector optics. Only the part of the beam that is scattered within a sample and returned 
to the surface is considered to be diffuse reflection. 

15 Few studies have investigated the potential of diffuse reflectance spectroscopy of 

tissue for detecting tissue cancerous changes. However, probes in the art used for diffuse 
reflectance measurements rely upon contact with the tissues to derive data that identify tissue 
that is suspected of being physiologically changed as a result of pre-cancerous or cancerous 
activity. Moreover, there has been no clear modelling approach that relates the changes in 

20 optical absorption and scattering coefficients to the physiological and morphological 
parameters to reflect early cancerous changes in real tissues. 

Despite the fact that employment of the contact probe geometry gives the most 
accurate diffuse reflectance measurements, the limitation of this kind of measurement is that 
it interferences with biopsy and other therapeutic procedures. Accordingly, the present 

25 invention uses a non-contact probe, which eliminates the need for a fiber probe through the 
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instrument channel, making clinical applications of this technology much more convenient. It 
also creates measurement geometry of broadbeam illumination and narrow spot detection, 
simplifying theoretical modelling of the measured spectra. 

FIG. 1 shows a system 12 for imaging and spectroscopy as is used in the present 
5 invention. The system has a light source 1, an endoscope 2 with a probe 3 that is adapted for 
insertion into the patient, a detection system including an image capture device 4 (such as a 
camera), a spectral attachment 10, and a spectrometer 5. The light source 1 provides 
illuminating radiation, preferably broadband light, via an optical fiber bundle 7 to an 
endoscope 2. The illumination fiber optic bundle 7 extends through the endoscope 2 and 

1 0 probe 3 to direct the illuminating radiation onto an investigated tissue 6. 

Light source 1 preferably is a Xenon arc lamp that provides both white light (or light) 
for white light imaging and reflectance spectral measurements, and a strong blue light (400- 
450 nm) with a weak near-infrared (NIR) light for fluorescence imaging and fluorescence 
spectral measurements. The NIR light is employed to form an NIR reflectance image used to 

15 normalize the green fluorescence image. (Another embodiment of the present invention 
employs a second excitation signal for a second fluorescence image used for normalization of 
the fluorescence image.) The light source 1 could also be a mercury lamp, a tungsten halogen 
lamp, a metal halide lamp, laser, or LED. Various filters may be added to select a given set 
of wavelengths. 

20 A processing unit 8 receives data from image-capture device 4 and spectrometer 5, 

and performs the calculations and processing as described herein. For example, the 
processing unit 8 will receive the diffuse reflectance spectra, as hereinafter described, and 
perform the analysis, classifying, and measuring functions described herein. The processing 
unit 8 is a computer or a microprocessor, preferably a personal computer. The processing 

25 unit 8 outputs its results to any output means desired by the user, such as a monitor, an LCD 
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screen, or a printer, or conveys the results to another computer for further analysis, or uses 
the results for its own internal calculations and analysis. 

The returning radiation from the tissue 6, which can be some combination of 
reflectance light, narrow-band emission light for fluorescence, other narrow bands for 
normalization or other types of light, is collected by various lenses and through the imaging 
bundle 9 is relayed to the detecting system for imaging by image detection device 4 and 
spectroscopy by spectrometer 5. 

Spectral measurements are performed using spectral attachment 10 mounted between 
the endoscope 2 and the detecting system. An optical fiber 11 carries a fraction of the 
returning radiation from a spot at an image plane to the spectrometer 5 for spectral analysis. 

The distal end of probe 3 is shown in more detail in FIG. 1 A. Endoscope probe 3, as 
shown in cross-section in FIG. 1A, typically contains one or more fiber optic illumination 
guides 21 to carry the interrogating radiation to the target object (such as the tissue 6 of FIG. 
1) and an imaging bundle 22 to cany returning radiation from the tissue 6. Probe 3 also 
contains an instrument channel 23 for biopsy or other surgical procedures, a water tube 24 for 
lavage of the target, and an air tube 25 for suction. In addition, the instrument channel 23 
may provide access for other medical procedures, such as optical computed tomography, 
Raman spectroscopy, confocal microscopy, endo-microscopy, laser or drug treatments, gene- 
therapy, injections, marking, implanting, or other medical techniques. A second imaging 
modality, such as fluorescence imaging, fluorescence spectroscopy, optical coherence 
tomography, Raman spectroscopy, confocal microscopy, or white-light reflectance imaging 
can be combined with the diffuse reflectance spectroscopy modality. 

In one embodiment, a light source 26 is placed near the distal end of the endoscope 
probe 3, as shown in FIG. 1 A. For example, by placing at least one LED and preferably at 
least two LEDs at the end of the endoscope, the fiber optics carrying the illumination or 
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excitation light can be eliminated. LEDs are lower cost, more reliable, longer lasting, lighter 
in weight, more compact and more efficient than lasers and lamps sources, allowing for better 
control of imaging and illumination. In addition, a miniature image-capture device 28 can be 
placed at the distal end of the endoscope 3, as shown in FIG. 1A. This configuration 
5 eliminates the need for a fiber optic bundle to channel the returning radiation to the image 
capture device. Instead, the miniature image-capture device 28 sends signals to a processor 
such a processing unit 8. This configuration provides an opportunity for increased resolution 
and improved imaging. In yet another embodiment the illumination sources, image detectors 
and other expensive optics may be disposed on a removable tip which can be changed from 
10 patient to patient, as described in detail in United States Patent Application No. 11/088,561, 
"Endoscopy Device with Removable Tip" the disclosure of which is incorporated herein by 
reference. 

FIG. IB shows the spectral attachment used in the endoscope system of the present 
invention. An endoscope 31 employing this spectral attachment 10 is described in detail in 

15 United States Patent 6,898,458, to Haishan et. aL, Methods and apparatus for fluorescence 
and reflectance imaging and spectroscopy and for contemporaneous measurements of 
electromagnetic radiation with multiple measuring devices^ the disclosure of which is 
incorporated herein by reference. This patent discloses various devices and configurations 
for simultaneous white-light and fluorescence imaging along with spectroscopy 

20 measurements. 

The light coming out from the illuminated tissue 6 of FIG. 1 is focused by lens 36 to 
form an interim image at the fiber-mirror assembly 32. The centre of the mirror is modified 
by drilling a hole in the center and an optical fiber 33 is inserted in the hole to take that 
fraction of the image to a spectrometer 35. The fiber position is seen in the image as a black 
25 spot indicating exactly where the spectral analysis will be carried out. The optical fiber 33 
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carries the reflectance signal from the spot of the image plane (point spectral measurement), 
which corresponds to an area of 1mm diameter at the tissue surface when the endoscope 
probe's tip is ten mm away from the tissue surface to the spectrometer 35 for spectral 
analysis. The physician can align the black spot with the area of interest and the spectral 
5 measurement along with the still image is saved into the computer memory. The video image 
and the spectrum (either in the WLR mode or in the FL mode) are simultaneously displayed 
on the computer monitor in live mode. Processing of returning radiation through multiple 
modalities, such as white-light imaging and fluorescence imaging, or imaging and 
spectrometry, is described in detail in the above-mentioned patent application and is 
1 0 incorporated herein by reference. 

A mirror 37 is placed in parallel to fiber-mirror assembly 32 to turn the light beam 
back to its original direction and then through lens 38 to a camera 34 for image acquisition. 
Mirrors 32 and 37 are at an angle of 45 degrees to the incoming light beam for illustration 
purposes only. 

15 The reflectance measurements performed by the system 12 of the present invention 

can be represented by an equivalent 1-D measurement geometry shown in FIG. 2. In such 
geometry, a continuous wave plane source irradiates the tissue and the reflectance is detected 
from a narrow spot on the tissue surface. The tissue, represented as a two-layer turbid media, 
is illuminated with a broadbeam S(z) of interrogating radiation through a non-contact 

20 perpendicular fiber optical bundle 7. The diameter of the illumination beam is approximately 
a two-cm spot on the tissue 6. When the tip of the endoscopy probe 3 is ten-mm from the 
tissue surface 6, for said diameter of the fiber bundle 7> we calculated that the returning 
radiation or reflectance signal is detected from a one-mm spot at the tissue surface 6. 

The method for diffuse reflectance measurement of the present invention will now be 

25 described. 
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The in vivo reflectance signal I m i(^) measured from the tissue can be described as 
following: 

5 where is the instrument spectral response, including source spectral features, fiber-bundle 
transmittance, and detector efficiency; ai is a constant related to the efficiency by which the 
tissue-surface specular reflection was collected by the probe; bi is a constant related to the 
efficiency by which diffuse reflectance from tissue is collected by the measuring probe; and 
{X) is the true tissue diffuse reflectance to be derived. 
10 In order to remove an instrument response, normalization should be performed using 

a standard reflectance disk with known reflectance. The process of removing of the 
instrument response is also shown diagrammatically in FIG. 3. 

The reflectance signal measured from the standard disc I m2 (A) , can be described as 
follows: 

15 I m2 (X) = a 2 I(X)+b 2 I(X)R s (2) 

where a 2 is a constant related to the efficiency by which the specular reflection is collected by 
the probe; b 2 is a constant related to the efficiency by which the diffuse reflectance is 

collected; and s is the reflectivity of the standard disc, which is a constant across the whole 
visible wavelength range and is very close to one. 
20 The in vivo reflectance signal Imi(^) measured from tissue is divided by the 

reflectance signal I m2 {X) measured from the reflectance standard disc to account for 
instrument spectral response I(X). Dividing Eq. 1 and Eq. 2 and rearranging the equation 
produces the following: 

KW = T^ = "o+boX<»W (3) 
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where, R nl (X) is the reflectance spectra measured by the apparatus 12 after removing the 
instrument response, R^iX)]^ the true tissue diffuse reflectance spectra, and a 0 and b 0 are 
additive offset and multiplicative factors respectively, which depend on the measurement 
conditions during each in vivo measurement performed. This includes the amount of specular 
5 reflection collected, the standard disc material used as reference, and the probe distance from 
the tissue during measurement. 

We have performed in vivo measurements on normal bronchial mucosa and both 
benign and malignant bronchial mucosa lesions on 22 patients and have obtained a total of 
100 spectra. A biopsy sample was obtained for each measurement to classify each measured 

10 tissue site into normal, benign or malignant. The pathology examination of biopsies revealed 
that 21 reflectance spectra were from normal tissue sites, 29 from benign lesions (26 
hyperplasia and three mild dysplasia), and 50 from malignant lesions (seven small cell lung 
cancer, three combined squamous cell carcinoma, 30 non-small cell lung cancer, ten 
adenocarcinoma). Our analysis was to develop algorithms to differentiate the spectra into 

15 two groups 

Group One: malignant lesions for tissue pathology conditions that were moderate 
dysplasia or worse; and 

Group Two: normal tissue/benign lesions for tissue pathology conditions that were 
better than moderate dysplasia. 
20 This binary classification is also in consistent with clinical practice that Group One 

lesions should be treated or monitored, while Group Two conditions could be left unattended. 
Also during routine clinical endoscopy examination, any suspected malignant lesions (Group 
One) should be biopsied, while Group Two conditions will not be biopsied. However, in this 
specially designed study, for each patient an extra biopsy was taken randomly from either a 
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normal-looking area or a suspected benign lesion so that we can assess the performance of 
the spectral diagnosis independent of the performance of the imaging diagnosis. 

We developed a method and device for tissue classification using quantitative 
information about tissue physiological and morphological changes obtained from the analysis 
of diffuse reflectance spectra. In order to achieve this we develop a forward model that relates 
tissue optical properties (absorption coefficient, scattering coefficient and scattering 
anisotropy) to the diffuse reflectance (computed) and than an Inversion algorithm to extract 
the quantitative information about tissue physiological and morphological properties from the 
tissue diffuse reflectance. 

Forward model 

The forward model is developed in the frame of light transport theory and discrete 
particle theory that relate computed diffuse reflectance Rc to specific tissue physiological and 
morphological parameters related to cancer changes. It is known that reflectance depends on 
the optical properties of the medium and the measurement conditions, such as distance from 
the probe, and angle in each measurement. For each tissue, light distribution in tissue can be 
described as a function of an absorption coefficient, a scattering coefficient, and a scattering 
anisotropy (direction of scatter). Therefore, we developed an absorption model with an 
optical absorption coefficient expressed in terms of the micro-vascular absorption-related 
parameters and a scattering model with a scattering coefficient expressed in terms of the 
tissue microstructure scatter-related parameters. 

FIG. 4 represents a block diagram of the forward model. The optical absorption 
coefficient (Inputl) is expressed in terms of blood volume content, oxygen saturation, in vitro 
lung optical properties (for example, the absorption coefficient of lung tissue measured in 
vitro with blood drained out) and oxy- and deoxy-hemoglobin absorption. The scattering 
coefficient (Input2) is expressed in terms of scattering volume fraction, scattering size 
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distribution function and tissue refractive index. All listed parameters are optical properties 
and are wavelength dependent. 

The present invention models a system with known optical parameters (absorption 
and scattering parameters) and calculates a computed value of the diffuse reflectance signal 
5 (Rc) in relation to those parameters. The light propagated in the tissue is modeled using the 
general diffusion approximation model used to analyze diffuse back-reflected light with 
estimated optical coefficients of tissue, which are correlated to morphological structure and 
chemical composition of the tissue. 

Theoretically, the forward model used to describe tissue reflectance spectra R C (X) at 
10 each wavelength, can be obtained using Fick's law: 



Rc(X) = 



= r l Vfcz,Jl)\^ 0 (4) 



where 0 is the light fluency spatial distribution, y is the diffuse flux, h is the incident power, 
and y is the diffusion constant, which depends on the tissue optical properties. The light 
fluency (j) was obtained from the general diffusion approximation model. The general 

15 diffusion equation is different from the standard diffusion approximation model in that it 
explicitly includes the collimated source in the radiance approximation and it uses the 5- 
Eddington approximation to model the single scattering phase function, and thus was 
expected to give better predictions of the visible light (470-700 nm) distribution in the 
superficial layer of lung, which was found to have low albedo value (i.e., \i a ~ p, s ). 

20 For a continuous wave plane source decaying exponentially in the z-direction, the 

general diffusion model is given by: 

VV(z)-/c,V(z) = -^(z) (5) 



20 



WO 2006/076810 PCT/CA2006/000080 

where <f> is the fluency rate; S(z) is the incident collimated source term, ju ir is the transport 
attenuation coefficient equivalent to [ju a + //y(l — where fi a and ju s are absorption and 
scattering coefficients respectively; ju* is the total attenuation coefficient and is equivalent to 
[ju a +//*]; //* is the reduced scattering coefficient which is equivalent to ju s (l- f) 9 where / 
is the fraction of light scattered forward in the 8-Eddington approximation to the scattering 
phase function; and g* denotes the degree of asymmetry in the diffuse portion of the 
scattering. The values of / and g* were related to the reduced single scattering anisotropy 
g from the matching of the second moment of the 5~Eddington phase function to the 

or 

Henyey-Greenstein phase function, and are equivalent to g and — — respectively. 

1+g 

Since we are interested in the diffuse reflectance that is more affected by the tissue 
mucosa superficial layer (approximately a 0.5-mm thickness), within which most of the early 
cancerous changes occur, we solved equation (5) for the two-layer tissue geometry (FIG.2) 
with the top layer thickness / — 0.5 mm. 

The light fluence distribution <f> as a function of tissue depth z (FIG.2a) is obtained 
using Monte Carlo simulation. The optical properties of lung tissue used in this simulation 
are obtained from Qu et ah described in detail in Optical Properties of Normal and 
Carcinomatous Bronchial Tissue, 33 Appl. Opt. 7397-405 (1994), the disclosure of which is 
incorporated herein by reference, and a four percent blood volume content is added into the 
tissue model. The fluence <f) becomes insignificant (reduced by factor e" 1 ) for depths after 0.8 
and 1 .6 mm for X = 470 nm and X = 700nm respectively. Therefore, the measured reflectance 
signal comes from a tissue volume starting from the surface and up to a depth of 0.8 to 
1 ,6mm depending on the wavelength of the light. 
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We solved Equation (5) in the z-direction for the (1-D) approximation model and the 
two-layer geometry, for layer 1 and layer 2, using the index mismatching boundary 
conditions at interface 221 (air-tissue interface) and the index matching boundary conditions 
at interface 222 (between the two tissue layers), as shown in FIG. 2. By substituting the 
solution of equation (5) into equation (4) we obtained an expression for the diffuse 
reflectance spectrum Rc{X) in terms of the absorption coefficient the scattering 

coefficient ju s , and the scattering anisotropy g. 

The absorption coefficient ju a 9 is modeled in terms of the blood contents and 
absorption coefficient of lung tissue measured in vitro with blood drained out. Two 
parameters that are used to describe the blood contents in tissue are the blood volume 
fraction p 9 and the blood oxygen saturations . The absorption properties of lung tissue in 
vivo can be described by the following equations: 

Ma W = Mblood WP + Mm vitro (1 ~ P\ ,~ 

Mbiood W = OMmoi + C 1 - a )Mm 
where jU HbQ2 and ju Hb are the absorption coefficients for oxy-hemoglobin and deoxy- 
hemoglobin respectively. The in vitro absorption coefficient vitro, is obtained from the in 
vitro lung tissue measurements made previously by Qu et ah described in detail in Optical 
Properties of Normal and Carcinomatous Bronchial Tissue,, 33 Appl. Opt. 7397-405 (1994), 
the disclosure of which is incorporated herein by reference. 

The scattering coefficient ju s and the scattering anisotropy g are modeled in terms of 
the microstructure scatterer volume fractions and size distribution. The tissue scattering 
model is developed using the fractal approach, assuming that the tissue microstructures 9 
refractive index variations can be approximated by a statistically equivalent volume of 
discrete ncdcro-scattering particles with a constant refractive index but different sizes. 
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For micro-scattering particles that are assumed to be spherical in shape, we can 
calculate the transport scattering coefficients for a bulk tissue by adding randomly the light 
waves scattered by each particle together. Thus, the transport scattering coefficient ju s and 
the scattering anisotropy g can be modeled using the following integral equations: 

5 M M=)\Q(^ri,X)]^dx (7) 

o U{X) 

Jfc(*,^)]^d* 

where Q(x) is the optical scattering cross section of individual particle with diameter x, 
refractive index n and wavelength X; v (x) is the volume of the scattering particle with 
diameter x and g(x) is the mean cosine of the scattering angles of single particle with diameter 
10 x. For spherical micro-particles, Q(x) and g(x) are calculated from Mie-theory using the Mie- 
scattering code described in details in C.F. Bohren and D.R. Huffman, Absorption and 
scattering OF LIGHT BY SMALL PARTICLES, the disclosure of which is incorporated herein by 
reference. 

The volume fraction distribution t](x) is assumed to follow a skewed logarithmic 
15 distribution: 



7j(x) = 8 C 0 x p exp 



(lnx-lnx w ) 



k 2 \ 



2aJ 



(9) 



J 



where 8 is the total volume fraction of all the scattering particles in tissue, fi is the size- 
distribution parameter (fractal dimension) which determines the shape of the volume-fraction 
size distribution and is related directly to the size of the scattering particles, x m and a m set the 
20 center and width of the distribution respectively, and C 0 is a normalizing factor obtained from 
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oo 

the condition S = [r[(x)(hc. The value of x m is assumed equal to the geometrical mean of 
0 

(0.05 jam) and (20 |im) which represent the limits of the scattering particles' range of 
diameters found typically in tissues. Thus x m = [(0.05)(20.0)] 1/2 = 1.0. The width parameter 
a m is assumed to be a constant of 2.0 to match with the fractal scaling range of tissues. 
Having x m and a m being set, the larger the value of p, the higher the contribution of the 
smaller size particles in the scattering particle size distribution function. 

The refractive index of the background surrounding media (nbkg) is assumed to be 
1.36. The refractive index of the scatterers inside the lung tissue is estimated based on the 
type of the tissue using the following relation: 

w= n bkg +ff(n r n s ) + {\-ff)(n n -n c \ 

(10) 

where n f is the refractive index of the collagen fibers which is equivalent to 1.47; n n is the 
refractive index of the nuclei which is equivalent to 1 .4; and n s and n c are the refractive index 
of the interstitial fluids and the intracellular fluids which are equivalent to 1.34 and 1.36 
respectively. The fibrous-tissue fraction ff value is assumed to be ten percent for the first 
layer (epithelial layer and part of the upper submucosal and is assumed to be 70 percent for 
the second layer, which is the lower submucosa and the cartilage layer. These assumptions 
result in refractive indexes of (nf=1.41) and (n 2 =lA5) for the first and the second layers 
respectively. 

Inversion Algorithm 

The Inversion (fitting) algorithm based on a Newton-type iteration scheme is 
developed to extract information about the tissue absorption and scattering parameters, taking 
into account measuring conditions (geometry correction parameters ao, bo), from the measured 
reflectance spectra. The Inversion consists of simulations with the different sets of 
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absorption and scattering parameters to determine the computed diffuse reflectance and 
determining tissue diffuse reflectance for a different set of geometry correction parameters 
(ao, bo) and than comparison between the simulations and the tissue diffuse reflectance. 
When we get the best match it allows us to determine the correct values of the absorption, 
scattering, and geometry parameters. 

The Inversion algorithm can be described through least-squares minimization 
function: 



where R m (/t.) is the reflectance measured at wavelength X t ; Rc(A l ) is the computed diffuse 

reflectance at wavelength X t according to Eq. (4), b Q is an intensity calibration factor to 

account for the instrument relative intensity measurements, and ao is an additive factor to 
account for the specular reflectance collected by the instrument probe during in vivo 
measurements. The following parameters are used as free-fitting variables during the 
inversion process: 

the blood volume fraction (p) assumed to be the same for both tissue layers, 

the blood oxygen saturation parameter (a) assumed to be the same for both tissue 

layers, 

the scattering volume fraction in the top layer (5\) and in the bottom layer (S 2 ) 5 

the size-distribution parameters (flj) and for both top and bottom layers 

respectively, and, 

the additive and multiplicative terms (geometry correction parameters) in Eq. 3 (ao ) 
and (b 0 ). 

Using the Gradient Base Search (Marquardt-type regularization scheme), we can 
obtain the updates of these parameters from the following system of equations: 




(11) 
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(C T C + "QAr = C T [R m - (a 0 + « )], 

(12) 

Where ^is the Jacobian matrix, At is the vector updates for the eight parameters (p, a, S u 
S 2 , fi i , P 2 , ao , bo), I is the identity matrix, and v is a scalar or diagonal matrix. The 
5 Jacobian matrix ^represents the sensitivity of the reflectance coefficients measured on the 
eight parameters and its elements are computed from the derivatives of Eq. (4) with respect to 
these eight parameters. The inclusion of a 0 and b 0 in the fitting are essential to account for 
the specular reflection component and the back-scattering probe collection efficiency, which 
vary for each measurement and depend, among others, on the probe-tissue distance and angle 
10 in each measurement. Thus, the true tissue diffuse reflectance R^W can then be extracted 
from the measured reflectance spectra R m (sl) measured by the apparatus, using the values of 
a 0 and b 0 that are obtained from the fitting procedure and substituting in equation (3). Of 
course, all other parameters, p, a 8\* S 2 , P i and ^2 are also derived by this Inversion 
algorithm. 

15 FIG. 5 illustrates the procedure disclosed with the present invention. The forward 

model and Inversion algorithm described above are used to derive the true diffuse reflectance 
spectra and to extract the cancer related physiological and morphological properties of the 
tissue under investigation. We ran a simulation using the inverse algorithm, with a set of 
known values of tissue optical parameters (p, a, S, ft) and compared the computed diffuse 

20 reflectance with the true tissue diffusion reflectance derived from the measured reflectance 
for known measuring conditions, (known geometrical correction parameters a 0 , bo). Once we 
get the best match, we had the actual values of tissue optical parameters related to diffuse 
reflectance. If there is no match, we ran another simulation with new set of parameters. We 
used a gradient base search to find new iteration parameters. 
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Statistical Analysis 

The actual parameters obtained using the Inversion algorithm are used for statistical 
analysis. The statistical analysis is performed to evaluate the differences in means between 
the two groups (benign and malignant), to determine which variables discriminate between 
5 the two groups, and finally to build classification functions for the evaluation of the 
developed model predictive classification. 

All fitted results obtained from the 100 spectral measurements are collected and saved 
in groups for statistical analysis. Because we have not been sure if the derived parameters 
follow normal distributions, the Kolmogorov-Smirnov two-sample test is chosen to evaluate 
10 the significance of the differences between the two groups (normal/benign tissue vs. 
malignant lesions) for each of the six parameters (p, a, S l:) S 2r> Pi) obtained from our fitting 

results. Discriminant function analysis (DFA) is then applied to the identified diagnostically 
significant parameters to build diagnostic algorithms for tissue classification. DFA determine 
the discrimination function line that maximized the variance in the data between groups while 

15 minimizing the variance between members of the same group. The performance of the 
diagnostic algorithms rendered by the DFA models for correctly predicting the tissue status 
{i.e. normal/benign vs. malignant) underlying each parameter set derived from the 
reflectance spectrum is estimated in an unbiased manner using the leave-one-out, cross- 
validation method on the whole data set. In this method, one case is removed from the data 

20 set and the DFA based algorithm is redeveloped and optimized using data of the remaining 
cases. The optimized algorithm is then used to classify the withheld spectrum. This process 
is repeated until all withheld cases (100 spectra/cases) are classified. The sensitivity and 
specificity are calculated from the results of the classification using the following 
expressions: 

25 Sensitivity = % (True Positives - False Negatives) /True Positives 

27 



WO 2006/076810 PCT/CA2006/000080 

Specificity = % (True Negatives - False Positives) /True Negatives 
The results of the statistical analysis performed to evaluate the differences in the 
means of the six measured parameters between benign and malignant group are summarized 
in Table 1 . 

5 Table 1 



Parameter 


Benign 


Malignant 


Significance (p) 


Mean 


Std. 
Dev. 


Mean 


Std. Dev. 


P 


0.032 


0.02 


0.065 


0.03 


0.001 


a 


0.9 


0.11 


0.78 


0.13 


0.022 


Si 


0.077 


0.057 


0.048 


0.046 


0.013 


Pi 


0.97 


0.15 


0.91 


0.12 


0.095 


8 2 


0.066 


0.048 


0.07 


0.032 


0.25 


P2 


0.94 


0.12 


0.92 


0.1 


0.65 



As shown in the table, the mean value of the blood volume fraction is higher for malignant 
lesions (0.065±0.03) compared to the benign lesions (0.032± 0.02). The mean value of the 
oxygen saturation parameter is reduced from 0.9 (for benign lesions) to 0.78 for malignant 

10 lesions. For the scattering related parameters, the mucosal layer show moderate to significant 
changes between normal/benign tissues and malignant lesions for the top layer with the mean 
values of Sj and J3j for the benign lesions to be 0.077 and 0.97 respectively, compared to 
0.048 and 0.91 for the malignant lesions. The scattering parameters (S 2 and Pi) for the 
bottom layer show minimal differences between the norma^enign tissue and malignant 

15 lesions. It should be noted that the larger the value of p, the higher the contribution of the 
smaller size particles in the scattering particle size distribution function. Thus, an increase in 
p value indicates a decrease in the scattering particle average size. The statistical analysis, 
using the Kolmogorov-Smirnov two-sample test, show that the malignant group has 
significant increase in the blood volume fraction, p (p =0.001 <0.05), significant decrease in 

20 the oxygen saturation parameter, a (p =0.022 < 0.05), and significant decrease in the mucosa 
layer scattering volume fraction, di (p =0.013 < 0.05) compared to the benign group. The 

28 



WO 2006/076810 PCT/CA2006/000080 

results also show moderate significant decrease in the size-distribution parameter of the 
mucosa layer (fii) in the malignant group compared to the benign group (p =0.095 <0.1). 

It should be noted that the significant increase in the blood volume fraction of the 
malignant lesions measured in our study agreed with the biological observations that tumors 
5 and cancerous tissues exhibit increased microvasculature and accordingly increased blood 
content. The significant decrease in the blood oxygenation in the malignant lesions is 
consistent in that hypoxia-related changes are occurring during cancerous development and 
could be related to the increase in tissue metabolism rate, the lower quality of the tumoral 
microcirculation, and to the high proliferation rate of the cancerous cells. The significance 

10 decrease in the scattering volume fraction found in the measured malignant lesions is 
consistent with the results obtained by Bard et dl. for the lung cancer lesions and the results 
obtained by Feld et at for the colon polyps. The explanation for such decrease in the 
scattering volume faction is still poorly understood due to the complex nature of tissue 
scattering process. However, this decrease may be related to the decrease in the 

15 mitochondrial content, which has been found to contribute most significantly to the light 
scattering in the backward (reflectance) directions or to the changes in the refractive index of 
the cytoplasm due to an increased protein and enzyme content. The size-distribution 
parameter (/?/) decrease means increased scatterer particle sizes on average for malignant 
tissues as compared to normal/benign tissues. This is consistent with the fact that cancerous 

20 cells have larger nuclei than normal and benign cells. 

The results obtained by analyzing the 100 reflectance spectra measured in vivo using 
the above-described model and curve fitting are shown in the FIGS. 6 through 10. An 
example of the reflectance spectra measured in vivo from two malignant lesions (two non- 
small cell lesions) and two benign sites (one normal lesion and one hyperplasia-diffused 

25 lesion) from the same patient are shown in FIG. 6a. As can be seen, the measured reflectance 
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spectra have large intensity differences, which are related to the variations in the specular 
reflectance signal and the distance between the endoscope tip and the tissue surface for 
different measurements. The accuracy of the model fitting to the measured reflectance 
demonstrated the validity of the developed method. The true tissue diffuse reflectance 
5 spectra, Rtm(X) is then derived by correcting the in vivo measured reflectance spectra R m (X) 
using the fitting results. The corrected true tissue diffuse reflectance spectra, Rtmft), are 
shown in FIG. 6(b). The specular reflection components have been successfully removed and 
the reflectance intensities fall between 0 and 1 (or 100 percent), while the original 
uncorrected spectra, R m (X)^ in FIG. 6(a) have quite arbitrary reflectance intensities between 0 

10 and 350 percent. The fitting results obtained from the analysis of the two benign and the two 
malignant spectra are summarized in FIG. 7. 

The average of the corrected reflectance spectra (Rtm) for both the normal/benign 
group and the malignant group are shown in FIG. 8. It shows that the average reflectance 
spectra of the normal/benign group have higher intensities in the measured wavelength range 

15 (470 - 700 nm) than the malignant group. This intensity differences are significantly larger 
for wavelengths above 600 nm. In addition the two hemoglobin absorption valleys around 
550 nm and 580 nm are larger and more obvious on the normal/benign group spectral curve 
than on the malignant group spectral curve. The average fitting parameters (p, a, 5i ,pi, 62, 
P2) and their standard deviations for the two groups are also shown in Table 1 . 

20 FIG. 9 shows the values of the bronchial mucosa layer (top layer) parameters (p, a, 

S l7 pi) obtained from the analysis of the 100 benign and malignant reflectance spectra 
measured. 

The results obtained from the DFA showed that the three parameters (p, a, S t ) are 
significant for the discrimination between the two groups. FIG. 10a shows the classification 
25 results based on measuring the blood volume fraction (p) and the scattering volume fraction 
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(dj) and (FIG. 10b) shows the classification results based on measuring the blood volume 
fraction and the tissue oxygen saturation parameters. As shown in the Figures, we can easily 
identify two domain spaces, with slight overlap, for benign and malignant groups. The DFA 
results with the leave-one-out, cross-validation method show that we could differentiate the 
5 measured lesions into normal/benign and malignant with sensitivity and specificity of 83 
percent and 81 percent respectively. Therefore, the method of the present invention, 
including illuminating tissue 6 with broadband interrogating radiation, collecting the 
returning radiation from the tissue 6 with a non-contact endoscope probe 3, measuring the 
diffuse reflectance spectra from the returned radiation and analysing measured diffuse 

10 reflectance spectra using the steps described above, results in classifying the tissue as benign 
or malignant with improved sensitivity and specificity. 

Please note that the method described above can also be used in conjunction with one 
or more of the imaging modalities described previously. 

While preferred embodiments of present invention are shown and described, it is 

15 envisioned that those skilled in the art may devise various modifications of the present 
invention without departing from the spirit and scope. 
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CLAIMS: 
We claim: 

1. A method of obtaining information about tissue physiology and morphology from 
diffuse reflectance spectra, comprising: 

illuminating a tissue with a broadbeam radiation to produce returning radiation; 
measuring a reflectance spectra of said returning radiation with a non-contact probe; 
determining a diffuse reflectance spectra from said measured reflectance spectra; 
analyzing said diffuse reflectance spectra for a two-layer tissue model by one- 
dimensional light transportation modelling; 

extracting at least one optical property of the tissue from said analyzed diffuse 
reflectance spectra; and 

deriving information about a physiology and a morphology of the tissue from said at 
least one said optical property. 

2. The method of claim 1 ? wherein said at least one optical property comprises at least 
one of an optical absorption coefficient, a scattering coefficient, and a scattering 
anisotropy. 

3. The method of claim 1, wherein said one-dimensional light transportation modelling 
comprises a forward model, an absorption model, a scattering model, and an inversion 
algorithm. 

4. The method of claim 3, wherein said forward model is used to model a system having 
known optical properties to calculate a computed value of said diffuse reflectance 
spectra. 

5. The method of claim 4, wherein said known optical properties are at least one of an 
optical absorption coefficient, a scattering coefficient, and a scattering anisotropy. 
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6. The method of claim 3, wherein said absorption model expresses an absorption 
coefficient in terms of blood contents and in vitro tissue optical parameters. 

7. The method of claim 6, wherein said absorption coefficient is at least one of an 
oxygen saturation and a blood volume fraction. 

8. The method of claim 3, wherein said scattering model expresses a scattering 
coefficient and a scattering anisotropy in terms of a scattering volume fraction and a 
size distribution parameter. 

9. The method of claim 8 ? wherein said scattering coefficient is at least one of a mucosa 
layer scattering volume fraction and a mucosa layer size-distribution parameter. 

10. The method of claim 3, wherein said inversion algorithm derives tissue physiological 
and morphological properties from said diffuse reflectance spectra. 

11. The method of claim 10, wherein said tissue physiological and morphological 
properties are at least one of an oxygen saturation, a blood volume fraction, a mucosa 
layer scattering volume fraction and a mucosa layer size-distribution parameter. 

12. The method of claim 3, wherein said inversion algorithm further derives values of 
geometry correction parameters. 

13. The method of claim 12, wherein said geometry correction parameters are used to 
determine said diffuse reflectance spectra from said measured reflectance spectra. 

14. The method of claim 12, wherein said geometry correction parameters account a 
specular reflectance collected by said non-contact probe and a variable collection 
efficiency of said non-contact probe. 

15. The method of claim 1, further comprising classifying the tissue as one of benign and 
malignant based on said optical property. 

16. The method of claim 15, wherein said classifying step further comprises comparing 
said at least one optical property to a data set of known pathology. 
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17. The method of claim 16, wherein said comparing step comprises using statistical 
analysis. 

18. The method of claim 1, further comprising at least one other modality for at least one 
of imaging and spectroscopy. 

19. The method of claim 18, wherein said at least one other modality is chosen from the 
group consisting of fluorescence imaging, fluorescence spectroscopy, optical 
coherence tomography, Raman spectroscopy, confocal microscopy, or white-light 
reflectance imaging. 

20. The method of claim 1, further comprising imaging the tissue with fluorescence 
imaging comprising at least a first excitation signal for a first fluorescence image and 
a second excitation signal for a second fluorescence normalization image. 

21. The method of claim 1, wherein said measuring step measures said reflectance spectra 
from substantially a point in an image plane of said returning radiation. 

22. An apparatus for obtaining information about tissue physiology and morphology from 
diffuse reflectance spectra, comprising: 

means for illuminating a tissue with a broadbeam radiation to produce returning 
radiation; 

a non-contact probe to measure said returning radiation; 

means for measuring a reflectance spectra of said returning radiation; 

means for detemiining a diffuse reflectance spectra from said measured reflectance 

spectra; 

means for analyzing said diffuse reflectance spectra for a two-layer tissue model by 
one-dimensional light transportation modelling; 

means for extracting at least one optical property of the tissue from said analyzed 
diffuse reflectance spectra; and 
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means for deriving information about a physiology and a morphology of the tissue 
from said at least one optical property. 

23. The apparatus of claim 22, wherein said at least one optical property comprises at 
least one of an optical absorption coefficient, a scattering coefficient, and a scattering 

5 anisotropy. 

24. The apparatus of claim 22, wherein said one-dimensional light transportation 
modelling comprises a forward model, an absorption model, a scattering model, and 



an inversion algorithm. 

25. The apparatus of claim 24, wherein said forward model is used to model a system 
10 having known optical properties to calculate a computed value of said diffuse 

reflectance spectra. 

26. The apparatus of claim 25, wherein said known optical properties are at least one of 
an optical absorption coefficient, a scattering coefficient, and a scattering anisotropy. 

27. The apparatus of claim 24, wherein said absorption model expresses an absorption 
15 coefficient in terms of blood contents and in vitro tissue optical properties. 

28. The apparatus of claim 27, wherein said absorption coefficient is at least one of an 
oxygen saturation and a blood volume fraction. 

29. The apparatus of claim 24, wherein said scattering model expresses a scattering 
coefficient and a scattering anisotropy in terms of a scattering volume fraction and a 

20 size distribution parameter. 

30. The apparatus of claim 29, wherein said scattering coefficient is at least one of a 
mucosa layer scattering volume fraction and a mucosa layer size-distribution 
parameter. 

31. The apparatus of claim 24, wherein said inversion algorithm derives tissue 
25 physiological and morphological properties from said diffuse reflectance spectra. 
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32. The apparatus of claim 31 ? wherein said tissue physiological and morphological 
properties are at least one of an oxygen saturation, a blood volume fraction, a mucosa 
layer scattering volume fraction and a mucosa layer size-distribution parameter. 

33. The apparatus of claim 24, wherein said inversion model further derives values of 
geometry correction parameters. 

34. The apparatus of claim 33, wherein said geometry correction parameters are used to 
determine said diffuse reflectance spectra from said measured reflectance spectra. 

35. The apparatus of claim 33, wherein said geometry correction parameters account a 
specular reflectance collected by said non-contact probe and a variable collection 
efficiency of said non-contact probe. 

36. The apparatus of claim 22, further comprising means for classifying the tissue as one 
of benign and malignant based on said tissue optical property. 

37. The apparatus of claim 36, wherein said classifying step further comprises comparing 
said at least one optical property to a data set of known pathology. 

38. The apparatus of claim 37, wherein said comparing step comprises using statistical 
analysis. 

39. The apparatus of claim 22, further comprising means for at least one other modality 
for at least one of imaging and spectroscopy. 

40. The apparatus of claim 39, wherein said at least one other modality is chosen from the 
group consisting of fluorescence imaging, fluorescence spectroscopy, optical 
coherence tomography, Raman spectroscopy, confocal microscopy, or white-light 
reflectance imaging. 

41. The apparatus of claim 22, further comprising means for imaging the tissue with 
fluorescence imaging comprising at least a first excitation signal for a first 
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fluorescence image and a second excitation signal for a second fluorescence 
normalization image. 

42. The apparatus of claim 22, wherein said means for measuring measures said 
reflectance spectra from substantially a point in an image plane of said returning 

5 radiation. 

43. A system for measuring quantitative information related to cancerous changes in a 
tissue, comprising: 

a non-contact probe; 

a light source producing a broadbeam interrogating radiation at a distal end of 
10 said probe to illuminate a tissue and to produce returning radiation; 

a detecting system coupled to capture said returning radiation; and 
a processing unit coupled to said detecting system, said processing unit 
measuring a reflectance spectra of said returning radiation, determining a 
diffuse reflectance spectra of said measured reflectance spectra and classifying 



15 the tissue as one of benign and malignant based on said diffuse reflectance 

spectra. 

44. The apparatus of claim 43 wherein said processing unit further measures at least one 
optical property of the tissue derived from said diffuse reflectance spectra. 

45. The apparatus of claim 44, wherein said at least one optical property comprises at 
20 least one of an optical absorption coefficient and a scattering coefficient. 

46. The apparatus of claim 45, wherein said at least one optical property comprises at 
least one of a blood volume fraction, an oxygenation saturation parameter, a mucosa 
layer scattering volume fraction, and a mucosa layer size-distribution parameter. 

47. The apparatus of claim 43, wherein said processing unit models a computed diffuse 
25 reflectance spectra for known optical properties using a forward model. 
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48. The apparatus of claim 43, wherein said processing unit extracts said at least one 
optical property of tissue from said diffuse reflectance spectra using an inversion 
algorithm. 

49. The apparatus of claim 48, wherein said inversion algorithm further derives values of 
geometry correction parameters. 

50. The apparatus of claim 49, wherein said geometry correction parameters are used to 
determine said diffuse reflectance spectra from said measured reflectance spectra. 

51. The apparatus of claim 49, wherein said geometry correction parameters account a 
specular reflectance collected by said non-contact probe and a variable collection 
efficiency of said non-contact probe. 

52. The apparatus of claim 43, wherein said processing unit further comprises means for 
comparing said at least one optical property to a data set of known pathology. 

53. The method of claim 52, wherein said means for comparing step uses statistical 
analysis. 

54. The apparatus of claim 43, wherein said light source comprises at least one of a 
Xenon arc lamp, a mercury lamp, a tungsten halogen lamp, a metal halide lamp, a 
laser, and a light-emitting diode. 

55. The apparatus of claim 43, wherein said light source further comprises wavelength 
filters. 

56. The apparatus of claim 43, wherein said light source is at a distal end of said probe. 

57. The apparatus of claim 43, wherein said probe is an endoscope. 

58. The apparatus of claim 43, wherein said probe further comprises a removable tip and 
wherein said light source is on said removable tip. 

59. The apparatus of claim 58, wherein said probe is an endoscope. 
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60. The apparatus of claim 43, further comprising at least one light guide coupling said 
light source to said probe. 

61 . The apparatus of claim 60, wherein said probe is an endoscope. 

62. The apparatus of claim 43, wherein said detecting system comprises at least a 
spectrometer. 

63. The apparatus of claim 43 , wherein said detecting system comprises an image capture 
device and a spectrometer. 

64. The apparatus of claim 63 , wherein said probe further comprises a removable tip and 
wherein said image capture device is on said removable tip. 

65. The apparatus of claim 43 further comprising an imaging bundle coupling said probe 
to said image capture device. 

66. The apparatus of claim 43, further comprising means for at least one other modality 
for at least one of imaging and spectroscopy. 

67. The apparatus of claim 66, wherein said means for at least one other modality is 
chosen from the group consisting of fluorescence imaging, fluorescence spectroscopy, 
optical coherence tomography, Raman spectroscopy, confocal microscopy, or white- 
light reflectance imaging. 

68. The apparatus of claim 43, further comprising means for imaging the tissue with 
fluorescence imaging comprising at least a first excitation signal for a first 
fluorescence image and a second excitation signal for a second fluorescence 
normalization image. 

69. The apparatus of claim 43, wherein said means for measuring measures said 
reflectance spectra from substantially a point in an image plane of said returning 
radiation. 
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Figure 2a: Light fluence distribution as a function of tissue depth , 
obtained using Monte Carlo simulation in a turbid media with 
optical properties typical to the lung tissue with 4% blood added 
( atA, = 470nm s and at ^=700 nm) 
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Figure 6: (a) Reflectance spectra measured (and fitted) from two normaVbenign tissue sites and two 
malignant lesions from the same patient. Nl - normal, N2 - hyperplasia-difnised, CI — non-small cell 
carcinoma, C2 - another non- small cell carcinoma. The symbols are experimental data points, while the 
curves are model fitting, (b) True tissue reflectance spectra derived from Eq. (3) using the correction 
coefficient do and bo obtained from curve fitting with our developed model. 
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Figure 7: Absorption and scattering related parameters obtained from fitting the two benign a bing the 
scatterer's size-distribution in both layers (fii and #?). nd two malignant spectra: (a) Blood volume fraction 
(p), (b) Tissue oxygen saturation parameter (of), (c) Scattering volume fraction in both layers (Sj and d 2 \ 
and (d) Size-distribution parameter descri 
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Figure 8: Average true reflectance spectrum of the 50 normal tissue/TDenign lesions versus that of the 50 
malignant lesions. 
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Figure 9: Scatter plots of the physiological and morphological parameters of the bronchial mucosa layer 
obtained from the reflectance spectral analysis: (a) Blood volume fraction p and (b) Oxygen saturation 
parameter a 9 (c) Scattering volume fraction^ , and (d) Size-distribution parameter/? 
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Figure 10: Binary plot of (a) blood volume fraction p versus the scattering volume 8± 9 (b) blood volume 
fraction p versus oxygen saturation parameter a. 
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